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Dark matter (DM) annihilations in the Galaxy may produce high energy neutrinos, 
which can be detected by the neutrino telescopes, for example IceCube, ANTARES 
and Super-Kamiokande. The neutrinos can also arise from hadronic interaction be- 
tween cosmic ray and atmosphere around the Earth, known as atmospheric neutrino. 
Current measurements on neutrino flux is consistent with theoretical prediction of 
atmospheric neutrino within the uncertainties. In this paper, by requiring that the 
DM annihilation neutrino flux is less than the current measurements, we obtain an 
upper bound on the cross section of dark matter annihilation (o~v). Compared with 
previous investigations, we improve the bound by including DM substructure con- 
tributions. In our paper, two kinds of substructure effects are scrutinized. One is 
the substructure average contribution over all directions. The other is point source 
effect by single massive sub-halo. We found that the former can improve the bound 
by several times, while the latter can improve the bound by 10 1 ~ 10 4 utilizing the 
excellent angular resolution of neutrino telescope IceCube. The exact improvement 
depends on the DM profile and the sub-halo concentration model. In some model, 
IceCube can achieve the sensitivity of (av) ~ 10~ 26 cm 3 s~ 1 . 

PACS numbers: 13.15. +g, 95.35. +d, 95.55.Vj, 98.62.Gq 
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I. INTRODUCTION 

Many astronomical observations indicate that most of the matter in our universe is dark 
(see e.g. Ref. [li]). The evidences come mainly from the gravitational effects of the dark 

nn 

matter (DM), such as the rotation curves of spiral galaxies [2j, |3j, the gravitational lensing 
4] and the dynamics of galaxy clusters 5|. The studies such as primordial nucleosynthesis 



6|] and cosmic microwave background (CMB) 



71 



show that the DM is mostly non-baryonic. 



Combining recent cosmological measurements, for example from the Wilkinson Microwave 
Anisotropy Probe (WMAP), one could deduce precisely the relic density of DM, namely 
&DMh 2 = 0.1143 ± 0.0034 8|. However, the nature of dark matter is still unclear. In the 
literature there is a "zoo" of particle candidates for DM 9[, among which the most popular 
candidate at present is the weakly interacting massive particle (WIMP) such as the lightest 
super symmetric particle (LSP), lightest Kaluza-Klein particle (LKP) etc. 

Search for WIMP in particle physics experiments is very important to pin down the prop- 
erties of the DM. Besides searching missing energy signals at accelerator-based experiments, 
there are usually two classes of methods to detect WIMP, namely direct and indirect ones. 
The former method detects WIMP by measuring the recoil of heavy nucleus in the detector 
and gives the most strong evidence for the existence of DM. The latter one detects the DM 
self-annihilation signals, which include neutrinos, photons, anti-matter particles and so on. 
Among them neutrinos are one of the most attractive signals. For the low energy neutrinos 
(say much less than 100 GeV), their interactions with matter are highly suppressed by a 
factor at least Q 2 /m 2 v with Q the typical energy scale of the interaction. The neutrinos are 
hardly energy loss and trajectory deflection during their propagation, therefore they may 
carry the information of the nature and distribution of the DM. However due to the same 
reason, it is hard to capture such kind of low energy and relatively low flux neutrinos. For 
the high energy neutrinos (say around 100 GeV or higher), the interactions among neutrinos 
and matter become much stronger. These neutrinos may keep the information of the DM, 
and it is relatively easy to observe them experimentally. 

One proposal of detecting the high energy neutrino signals is to explore the locations 
close to us such as the center of the Sun or the Earth to get enough neutrino flux. The DM 
particles are grav itational trapped in the center of the Sun or Earth and produce neutrinos by 



annihilation 



lOt ] . If the annihilation and capture processes are in equilibrium, the neutrino 
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flux are mainly determined by the cross section of the DM and nuclei. Another proposal 
is to detect the neutrino signals from DM annihilation in the Milky Way (MW). Though 
the sources in the MW are farther than the Sun, it is natural to expect that the larger 
amount of DM can compensate the distance. Moreover the neutrino flux depends on the 
DM annihilation rate and number density square, therefore the regions with high density 
in the MW, such as the Galactic Center (GC) or sub-halos, should be potentially excellent 
observational targets. 

The GC is conventionally thought to be source-rich astrophysical laboratory and has at- 
tracted many attention of astronomers. It is also true for the DM indirect searches. Due to 
the weak interaction of DM particles, the DM density at the GC is highly accumulated as 
shown by detail simulations, which makes the GC a bright source of DM annihilation. How- 
ever, the complicated astrophysical environment and various kinds of astrophysical sources 
make the GC a high background site. In addition, the overlapping with bayonic matter 
objects (e.g., the central massive black hole ll|) may also affect the DM distribution and 
increase the uncertainties. It should be emphasized that the galactic sub-halos may be good 
candidates as DM sources. Since the self-annihilation of DM is square-dependent on the 
number density, the clump of substructure is expected to effectively enhance the annihi- 
lation signal and plays a role of the so-called "boost factor" 12, Q]- Furthermore as 



indicated by simulations, the spatial distribution of DM sub-halos is tend to be spherical 
symmetric in the MW halo, which may locate at a low-background site and effectively avoid 
the source confusion in the galactic plane. The effects of DM sub-halos on the flux of induced 
neutrino is what we try to investigate in this work. 

The neutrinos detected by high-energy neutrino telescopes such as Super-Kamikande 



AMANDA 



161 ] etc. are thought to be mainly from atmospheric neutrinos. Here they 



originate from the decay of hadrons which are produced by the strong interactions of cosmic 
rays with atmosphere. Experimentally no obvious excess has yet been observed. Then the 
measurements on neutrino flux can be utilized to set bounds on the DM annihilation cross 
section. As neutrinos are the most difficult to detect in the DM annihilating final states, 
the authors of Ref. and 18] assumed that the DM annihilate solely into neutrinos. 
They calculated the extra-galactic and the galactic neutrino fluxes, compared them with 
atmospheric neutrino flux, and set an upper bound on the DM annihilation cross section 
(crv). However, the DM may annihilate into final states other than neutrinos. In most of 



i 



the DM models there are several annihilation channels. Moreover high energy neutrinos 
from the DM annihilation will lead to gauge bosons bremsstrahlung [igL I20I even in the 



standard mode 
also inevitable 



(SM). The electromagnetic final states through higher-order corrections are 



2l| . Therefore the assumption that DM annihilate into only neutrinos gives 



the most conservative bound on DM annihilation cross section. 

In this paper we calculated the neutrino flux from the DM annihilations in the MW 
including the contributions from sub-halos by assuming that the DM annihilate into neutri- 
nos only [n], Q|- By comparing the predicted flux with the available atmospheric neutrino 
measurements, we set a very strict constraint on the DM total annihilation cross section. 
Compared to the previous studies, in this work we utilize the angular resolution of the neu- 
trino telescope to derive the stricter constraints. Here the massive sub-halos can be treated 
as the point-like sources. Based on our analysis we may observe the hig h energy neutrino 



22 



and IceCube 



23] 



flux provided that precise angular resolution data from ANTARES 
is available. On the other hand if no excess flux out of atmospheric neutrino is observed, 
an improved upper-bound of the annihilation cross section and/or the exclusion of certain 
sub-halo models can be obtained. 

This paper is organized as following. In Sec. II, we describe the sub-halo models according 
to the N-body simulation results. In Sec. Ill, we give the constraints on the dark matter 
annihilation cross section. The conclusions and discussions are given in Sec. IV. 

II. GALACTIC DM DISTRIBUTION AND SUBSTRUCTURE 



The current knowledge of the DM spatial distribution is mostly from the N-body simula- 
tion. Navarro et al. 24| firstly proposed a universal DM profile (referred as "NFW"). Based 



on their simulation in a wide range of halo mass, the density of DM can be written as 

{r/r s )[l + r/r s \ 2 

where and r s are the density scale and radius parameters for a specific DM halo. Moore 



et al. 



251 ] gave another profile with a more cusped inner slope compared with NFW as 

Ps 



p(r) 



{r/r s y- 



In addition, cored profile as 



(l+r/r.Vi 



(r/rs) 2 ] 
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(l/r s ) L5 ] 



(2) 



27| or cuspy profile with different inner 



slope from NFW and Moore (e.g., [28() were also proposed in the literature. Reed et al. 
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showed that even the inner slope steepens with the decrement of the halo mass, instead of a 
universal one 



29]. All these profiles give the same behaviors ~ r~ 3 at large radii, but show 
discrepancies in the inner region of the halo. Precise determination of the DM profile needs 
simulation with higher resolution and better understanding of the DM properties such as the 
interaction with baryonic matter. In this work we will adopt both NFW and Moore profiles 
for the discussion. It should be noted that the central density for NFW or Moore profile is 
divergent. In order to handle the singularity, a characteristic radius r c is introduced within 
which the DM density is kept a constant p max due to the balance between the annihilation 



rate and the in-falling rate of DM 
31]. 



30|. Typically we have p max = 10 18 ~ 10 19 M kpc" 



A. Determination of the profile parameters r s and p s 



Following 



321 ]. we use the virial mass M v of the halo to determine the parameters p s and 



r s . For a DM halo with specified mass M v , the virial radius r v is defined as 

(3) 



(4tt/3)A Pc 

with the density amplifying factor over the background A m 200 and the critical density of 
the universe p c = 139 M kpc~ 3 . The concentration parameter c v (CP) is defined as 

c v = (4) 
r-2 

d(r 2 p) 

where r„ 2 refers to the radius at which K dr ' | r=r _ 2 = 0. It is shown that for NFW profile 
r_ 2 = r™' w and for Moore profile r_ 2 = 0.63r™ oore , so we have 

nfw _ T _v_ 

(M v ) 

„moore _ 

(M v ) 

c v {M v y s 0.63 c v (M v y 1 ' 

Then p s can be derived just by requiring J p(r)dV = M v . We can see that the profile of the 
DM halo is fully determined provided that the c v — M v relation is specified. 

Generally the c v — M v relation is fitted from the numerical simulation. Here we will use 



two toy models of Eke et al. ( 33J, denoted by ENS01) and Bullock et al. ( 32j], denoted 
by B01), within which the DM halo forms based on the cosmo logical background density 
field. The CP predicted in these models increases with the decrement of the halo mass. 
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Such behavior is understandable in the frame of hierarchy structure formation, i.e. smaller 
halo forms earlier when the universe is denser than today. This behavior is confirmed at 
the cluster scale 



34 



351 ] . However, other studies showed agreement or disagreement with 



these two models, which indicate that we may not achieve the final understanding of this 
topic at present (see the discussion of Ref. 3l|). For the current work, these models are 
regarded as reference ones. We use the fitted polynomial form of the simulation at z = 



and extrapolate to low masses 311 ]: 



ln(c„) = ^2 C i 



i=0 



In 



Mv_ 

Mr 







(6) 



with M the mass of the Sun and 



C 



ENS01 



{3.14, -0.018, -4.06 x 10" 4 , 0, 0} 



and 



C 



B01 



{4.34, -0.0384, -3.91 x 10~ 4 , -2.2 x 10~ 6 , -5.5 x 10" 7 }. 



(7) 



(8) 



Fig. [T] shows c v as a function of the halo mass M v (see also Fig. 1 of Ref. 36j). For the 
very low mass, it is shown that c v becomes flat due to the fact that small objects tend to 
collapse at the same epoch. It should be noted that ENS01 and B01 models are for distinct 
halos in the universe. For the sub-halos within a host halo which is denser than the universe 
background, as we will discuss in this work, it is expected to be more concentrated than the 



distinct halos. In Ref. 



321 ] the simulation showed the sub-halos within a host halo indeed 



have larger c v than the distinct ones with the same mass. The simulation shows the c v of 
subhalo has c v ~ M~ a3 , which is steeper than the distinct halo (c v ~ M~ a13 ). In Fig. [I] we 
also show the extrapolated results of c v for sub-halos with mass 10 6 ~ 10 12 M Q . The other 
way to deal with the sub-halo is to multiply the result for distinct one by an empirical factor 
2 in Ref. 



le.g-, 



361 ] . in the following we denote this model as S01 x 2). 



B. The MW halo and substructure 



The mass of the MW DM halo is about 1 ~ 2 x 10 12 M determined from the rotation 
curve or kinematics of tracer populations such as the stars, satellite galaxies and globular 



clusters 



37 



38 



39j. A recent work by Xue et al. showed that M v = (0.93 ± 0.25) x 10 12 M G 



o 
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SUN 



FIG. 1: Concentration parameter (CP) c v as a function of halo mass M v at epoch z = 0. 



through an analysis of kinematics of a large sample of SDSS halo stars 40j. Here we adopt 
M v = 10 12 M as the mass of the total MW halo. As shown below, about 10% ~ 20% of 
the mass will form substructures, so the mass of the smooth halo is 0.8 ~ 0.9 x 10 12 M . A 
NFW profile is adopted for the smooth halo. Using the B01 model, we find that 1 r v pa 205 
kpc, c v ~ 13.6 and p Q ~ 0.34 GeV cm -3 . This configuration of the smooth halo is fixed in 
the following discussion since we will focus on the substructure in the present work. 



The surviva 



41 



42 



43 



44 



of substructure in galactic halo was revealed by many simulation groups 



45 



46 



471 ] . A recent simulation conduced by Diemand et al. showed 



that the self-bound substructure could even be as light as the Earth, with a huge number 
reaching ~ 10 15 28j. Simulations give the number density of sub-halos an isothermal spatial 
distribution and a power-law mass function as 

dA „ ( M sub 



dM sub ■ Anr 2 dr 



An 



Mi 



host 



1 + — 

r H 



2 ■ 



(9) 



1 Here 10 12 M Q is used to calculate r v and c v . while the density at the solar location is scaled to match 
that of the mass of the smooth component is 0.85 x 10 12 M©. 
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where M su & and Mh ost are the masses of sub-halo and host halo (10 12 M Q here), r# is the 
core radius which usually is a fraction of the virial radius of the host halo, a is the slope 
of the mass function and Nn is the normalization factor. For a galactic host halo, th was 



bund to be about 0.14r„ 



43 



44 



49 



50 



481 ] . The slope a lies between 1.7 and 2.1 in various works 



511 ] . Here we adopt a = 1.9 as in Ref. [12| 



14, 13611 - The mass function 



of Eq. Q is thought to be held in the mass range from the minimal sub-halo with mass 
~ 10 -6 M Q which is close to the free-streaming mass 28|, |52|, |53j, |54j], to the maximum one 
about 0.01 Mh ost [jjj]. The normalization is determined by setting the number of sub-halos 
with mass larger than 10 8 M is 100 3jJ. Finally in the inner region of the host halo, strong 
tidal force tends to destroy the sub-halo and the survival number should be cut down. We 



employ the "tidal approximation" as in Ref. 
mass fraction of substructure is about 14%. 



121 ]. Under this configuration, we find that the 



C. Astrophysical factor of the DM annihilation 

The annihilation signal of DM particles relies on two factors: the particle physical factor 
W(E) (energy dependent) depending on the particle property of DM, and the astrophysical 
factor J(ip) (spatial dependent) depending on the spatial distribution of DM. The neutrino 
flux observed on the Earth (applicable also for 7) can be written as 



9 1 (av) diV 1 



(p(E,iP) = C x W{E) x J{rj}) 

1 (av)dl 

4^2^f dE " p%R & J LOS ' 

where p Q = 0.34 GeV cm -3 is the local DM density and R Q = 8.5 kpc is the distance of 
the Sun from the GC, ip is defined as the angle between the observational direction and the 
GC direction relative to the observer, (av) is the average value of annihilation cross section 
times relative velocity, m x is the mass of DM particle, dN/dE is the production spectrum 
of v per annihilation. The integral path in Eq. (1101) is along the line-of-sight (LOS). To 
account for the contribution of substructures, we just need to replace p 2 in Eq. (fTOj) by 
(P 2 ) = pLooth + (PL) ™th Q 

dN ( r 2 



^ = / dM ^ dM^.4^dr U PLbdV 1 ' (U) 



s ■(.; b 



The average astrophysical factor within a solid angle AQ (e.g., the angular resolution of 
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the detector) is defined as 

JAn = ^y JWdfi. (12) 

where J(il>) is defined in Eq. ([TO]) . In Fig. [2] we show J(ip) as a function of ^ and Jah as a 
function of smooth angle AO, which is the half angle of the cone centered at the direction 
of the GC. Here AQ = 2ir(\ — cos AO). From the figures we can see that the enhancement 
on J(ip) by substructures is mainly at large angle, i.e., the direction deviation from the GC. 
The contribution from sub-halos with Moore profile is about 8 times larger than that of 
NFW profile (slightly different between various concentration models). For ENS01 model 
the enhancement upon the smooth component is very weak, while for the best combination 
B01x2+Moore the boost factor can be as large as ~ 25 at the anti-GC direction. However, 
the enhancement of Jaq is not remarkable. For the halo-average (AO = 180°) the best case 
gives ~ 4 times boost as shown in the right bottom figure in Fig. [2j 

D. DM substructure as point-like source 

In the previous section the contribution from substructures is averaged over the whole 
MW (see Eq. ( TTTT) ) . It should be noted that for the few massive sub-halos this average 
is unreasonable. In this case, the massive sub-halo should be more appropriately treated 
as point-like source and may be identified by high angular resolution detector such as the 
IceCube. This feature can also be utilized to suppress the background. 

In our work, a Monte-Carlo method is adopted to produce sub-halos with mass larger 
than 10 6 M according to the distribution function Eq. Qj. For this distribution function, 
we find N(> 10 6 M ) « 6400. In our numerical simulation, totally 50 MWs are generated, 
and the "tidal approximation" is adopted. For each sub-halo, we calculate the value of J{ip) 
and average within the cone with half angle equal to 1°. Then we count the cumulative 
number of sub-halos with astrophysical factor larger than a specified value, as shown in Fig. 
EJ The concentration models for ENS01, B01 and BOl-subhalo (see Fig. [lj are adopted. We 
can see that the ENS01 model gives the smallest astrophysical factor, while the result from 
BOl-subhalo model is larger by orders of magnitude than other models. From the figure 
we can see that the uncertainties of density distribution in sub-halos are very large. Given 
the particle factor of DM annihilation, we can get the number of sub-halos with flux higher 
than a specified value, such as the sensitivity of the detector. On the contrary, if no source 
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FIG. 2: J(ip) as a function of ip (left column) and Jah as a function of smooth angle AG — half 
angle of the smooth cone centered at the direction of GC (right column). In each figure the 
smooth (solid), smooth+sub(NFW) (dashed) and smooth+sub(Moore) (dotted) are shown. For 
each column, from top to bottom, the figures correspond to the concentration models as ENS01, 
B01 and B01 multiplying by a factor of 2, respectively. 
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is seen, it puts a constraint on the DM annihilation cross section 55] 
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FIG. 3: Cumulative number of sub-halos whose astrophysical factor > J±o as a function of J\c 
The dashed horizon line corresponds to the case that only one subhalo will be observed. 



III. CONSTRAINT ON DARK MATTER ANNIHILATION CROSS SECTION 
FROM MEASUREMENTS OF ATMOSPHERIC NEUTRINO FLUX 



The atmospheric neutrino comes from hadronic interaction between cosmic ray and atmo- 
sphere around the Earth. The observed neutrino flux is consistent with theoretical prediction 
within the uncertainties 0, l^]. The neutrinos from the DM annihilation should not be 
larger than the measured flux of the neutrino. Thus this requirement gives bounds on the 
DM annihilation cross section provided that the astrophysical factor, as discussed in the 
above section, is known. Actually in the literature the constraints on the neutrino fluxes 
from cosmic diffuse processes and the whole MW halo have been investigated in Ref. 13, Il8 |. 
In this work, influence on neutrino flux by DM substructures will be scrutinized. 

We assume that DM particles annihilate into pairs of neutrinos following Ref. 13]. The 
spectrum of the neutrinos per flavor is a monochromatic line with dN u /dE v = \8{E V — m 



X)i 



12 



where m x is the mass of the DM particle. All neutrino flavors are equally populated and 
the neutrino and anti-neutrino are added together. We require that the average neutrino 
flux from DM annihilation should not exceed the average atmospheric neutrino flux, which 



is taken from 



58]. 



We first consider the neutrino flux averaged over all directions with smooth and sub-halo 
contributions. As shown by the left column in Fig. [2J the sub-halo contribution to the 
neutrino flux is insensitive to the direction [3], while the smooth contribution can vary by 
several orders of magnitude for GC and anti-GC directions. The sub-halo contribution in 
anti-GC direction is larger than that of the smooth one, but smaller in GC direction. In other 
words, the neutrinos from anti-GC direction are mainly from sub-halo DM annihilation. 

Second, we investigate the effects on the neutrino flux from DM sub-halo contribution. 
Massive sub-halo is more appropriately treated as point-like source, as discussed above. We 
take smaller angle cone to calculate the neutrino flux from such single sub-halo and the 
atmospheric neutrino separately. The average neutrino flux from sub-halo in this cone is 
required not to exceed the atmospheric neutrino flux. Combined with the input astrophysical 
factor, we can set bounds on the DM annihilation cross section. This approach is obviously 
based on the excellent angular resolution of the neutrino telescope. The angular resolution 



of Super-Kamiokande is 69(E) ~ 30° x x/GeV/E |57[. For E > lOOGeV neutrino, the 
angular resolution can reach about 3° for Super-Kamiokande, while 1° for IceCube 23]. 
In our evaluation we adopt a 1° half-angle cone for IceCube. Note that for IceCube, the 
candidate source should be in the northern sky, and its threshold energy is 50GeV 



In order to obtain the constraint on the DM annihilation cross section, we must input the 
astrophysical factor. We use average neutrino flux in a cone </>ao by simply replacing J(i/j) 
with Jaq in Eq. ( flOl) . Since the DM induced neutrino is sharply peaked, we average neutrino 
flux within the energy bin around E v = m x to compare with the atmospheric neutrinos. 
And the energy bin width is taken as Alog 10 £' = 0.3 within the energy resolution limits 



of the neutrino detectors Super-Kamiokande [15|, |57|and IceCube [23j. The Jaq used in 
sub-halo point source is taken from Fig. [3] when cumulative number equals to one. Based 
on the astrophysical factor input this way, the derived annihilation cross section gives the 
sensitivity that IceCube may find at least one sub-halo. Jaq for different sub-halo profiles 
and concentration models are given in Table. [B 

From Tabled we can see that J^n is 3.0 for the smooth case. The sub-halo contribution 
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Halo average 


Point-like 




smooth smooth + subhalo 


ENS01 B01 BOl-subhalo 


NFW 


3.0 4.1 


5.3 28.9 12763.7 


Moore 


12.1 


43.3 276.3 82110.2 



TABLE I: Jac for different sub-halo profiles and concentration models. The values listed here 
correspond to the cumulative number equal to 1 as shown in Fig. [3j The smooth contribution to 
Jaq is fixed to be NFW profile. 'Halo average' means averaging contributions over the whole MW 
halo. For the sub-halo contributions we choose the concentration model of B01x2, as shown in 
the bottom-right figure in Fig. [2j ENS01, B01 and BOlsubhalo represent different concentration 
models when considering DM substructure as point-like sources. For point-like sub-halos, Jah is 
averaged in a 1° half-angle cone around the center of the sub-halo. 

with Moore profile is 9.1, while only 1.1 for NFW profile. If one averages the contributions 
over the whole sky, the sub-halo contribution can enhance the neutrino flux slightly. As 
shown in Fig. the enhancement of sub-halo is large at the anti-GC direction which can 
reach about 100 ~ 1000, but the enhancement at the GC direction is small. Note that the 
flux from the GC is much larger than that from the anti-GC, it is natural to expect the 
averaging sub-halo contribution should not be significant. Therefore the point-like sub-halo 
should be naturally more important, as shown in Table [H For the different concentration 
model ENS01, B01 and BOlsubhalo, Jac is much larger than the average contribution, 
though the uncertainty is very large. The numbers in different concentration models can 
be understood from Fig. CD, which indicates that the large c v corresponds to the large Jaq. 
Virial radius r v denned in Eq. H] is the boundary that gravitational force can sustain itself, 
and it is only related to virial mass M v but not the mass distribution. The density profiles 
with NFW or Moore distribution behave both like r -3 at large radii and r _1 or r~ L5 at 
small radii separately. The radius r_ 2 is a measure of the density profile p(r), namely the 
dark matter distributes mainly within the region r < r_ 2 . Thus larger c v means more dark 
matter concentrating in the sub-halo center, which results in a larger Jaq- 

Once the astrophysical factor is known, bounds on the dark matter annihilation cross 
section (av) can be derived by requiring the DM induced neutrino flux less than the measured 
ones. In Fig. [4] and Fig. [5j the upper bounds on (av) are plotted for two DM profiles. Also 
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shown in the figures is (av) for the natural scale of dark matter. The upper solid line 
correspond to the constraints from 180° half-angle cone averaging contributions. The other 
lines represent the sub-halo point source with 1° resolution from the IceCube. The dashed, 
dash-dotted and short-dashed lines are plotted with different concentration models. The 
constraints from the whole halo average are slightly improved compared with the smooth 
case m Refs. Q,Q- For the sub-halo as the point source, the constraints are significantly 
improved. Moreover the upper bounds of (av) are about 10~ 23 ~ 10~ 24 cm 3 s _1 for B01 
model, and can even reach ~ 10~ 26 cm 3 s _1 for BOl-subhalo model. Such bound is even 
lower than (av) for the natural scale which can induce the correct relic density of DM. 
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FIG. 4: The upper bounds on the DM annihilation cross section (av) as a function of dark matter 
mass with NFW subhalo profile. The upper solid lines represent the constraints from 180° half- 
angle cone average. The other lines represent constraints from sub-halo point source with 1° 
resolution from the IceCube. The dashed, dash-dotted and short-dashed lines are plotted for 
different concentration models. The lower solid line corresponds to (av) ~ 3 x 10 _26 cm 3 s~ 1 for 
the natural scale to produce the correct thermal relic density. 
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Subhalo with Moore distribution 
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FIG. 5: Same with Fig. S]but with Moore subhalo profile. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, by requiring the DM induced neutrino flux less than the measured ones, 
we give the improved upper bounds on the DM annihilation cross section (av ) with the DM 
substructure effects included. Here we assume the DM particles annihilate into neutrinos 
solely following the previous works. The observed neutrino flux depends on the particle 
physical and astrophysical factors. Thus we first investigate the astrophysical factor. Several 
different DM profiles and sub-halo concentration models are adopted based on the numerical 
simulations. Our studies show that at the anti-GC direction, the enhancement factor of sub- 
halos for B01x2 model is about 3 and 25 for NFW and Moore profiles respectively. While 
the whole-sky average (halo average, with cone half-angle 180°) does not have prominent 
enhancement. The best case of our adopted models is only ~ 4 times larger than the smooth 
ones (see bottom-right figure in Fig. [2]). If we take the 30° angular average (halo angular 
in Ref. [3]), there is almost no enhancement, as can be seen in Fig. [2j This is because the 
enhancement from sub-halos is spatially dependent on the MW halo, instead of a universal 
one On one hand, the smooth component increases more rapidly and dominates the 
annihilation flux near the GC; on the other hand, the tidal disruption on sub-halos is most 
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effective close to the GC. Thus the effect of substructures is not significant near the GC. In 
the cases of halo average and halo angular, the GC contributions are included and play 
a dominant role in the total flux, therefore no remarkable enhancement from sub-halos is 
found. 

In this paper we emphasize the important role of the massive sub-halos (e.g., M su f, > 
1O 6 M ). Since the number of massive sub-halo is small, it should be regarded as the point- 
like source. For the massive sub-halos, the high angular resolution of neutrino detector can 
be utilized to suppress the atmospheric neutrino background, thus the constraints on (crv) 
are expected to be improved. The angular resolution ~ 1° for energy greater than 50 GeV 
of the forthcoming experiment IceCube is employed 2 . The neutrino signal flux in a cone 
with half-angle 1° is calculated and compared with the atmospheric background in the same 
cone. We found that the constraints on (crv) are indeed improved significantly. Note that 
the constraints are model-dependent. For the moderate case B01+NFW, we find the upper 
bound of (av) is about 10 -23 cm 3 s _1 . While for the concentration model BOl-subhalo, the 
bound can reach 10~ 26 cm 3 s" 1 , which is even lower than (av) for the natural scale which 
can induce the correct relic density of DM. 

It should be noted that DM can annihilate into final states other than neutrinos. Thus 
the assumption that DM annihilates into only neutrinos gives the most conservative bound 
on the DM annihilation cross section and is independent of the particle properties of DM 
particle. 

Neutrinos are thought to be an important complementary particles for DM indirect 
searches besides 7-rays and charged anti-particles. It is shown that the detectability of 7- 
ray s from DM sub-halos on GLAST is optimistic |s9] . The effects of subhalos on positrons 
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6l| and antiprotons 31] are also investigated. The combination and cross check of dif- 



ferent kinds of signals will be very crucial to identify the DM sources and investigate the 
properties of DM particles. 



2 For Super-Kamiokande, the resolution angle is ~ 3° for E > 100 GeV and the background is 9 times 
greater. 
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